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Abstract 

We explain how spectrally stable vortices of the Nonlinear Schrodinger Equation in the plane 
can be orbitally unstable. This relates to the nonlinear Fermi golden rule, a mechanism which 
exploits the nonlinear interaction between discrete and continuous modes of the NLS. 


1 Introduction 

In this paper we consider the nonlinear Schrodinger equation (NLS): 

lut = —Am + Vu + /3(|mP)u, u(0, x) = uo{x), {t, a:) € M x (1-1) 

with V a real valued Schwartz function. We are interested in bound states, which are solutions of 
(1.1) of the form u{t,x) = e“^‘(/)tj(x) with w S R. When (j)^^ is real valued and of fixed sign then we 
call a ground state. In all other cases we call it an excited state. In the d = 2 case on which 

we focus in this paper, and with V{x) = V(|x|), an important class of excited states, which we call 
vortices, involves solutions of the form 

u{t,x) = e^^*<j)^{x) with (ji^iix) = with —>• R, (1.2) 

with <j)u){x) smooth and rapidly decreasing to 0 at infinity and m > 1. In the sequel, we will always 
assume that there is a family of bound states for ui in some open interval O C R+ = (0, oo) (see 
(H4) in section 2). We will study the following classical notion of stability, [31, 81]. 

Definition 1.1 (Orbital stability). A bound state of (1.1) is orbitally stable if 

Ve>0, 3(5>0 s.t. \\4)uj — wollui < <5 sup inf — u{t)\\Hi < e 

t>0 

where u is the solution of (1.1) with m(0) = uq. 

The orbital stability of bound states has been extensively studied, mainly using two tools: 
Lyapunov functions; linearized operators. 

It is well known that (1.1) conserves the energy 

E{u):=]-[ \yu\^ + V\u\^ + B{\u\^)dx (1.3) 

2 JR2 
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where B{0) = 0 and B'{s) = /3(s), and the mass 

/ d 

\u\^dx (we will set q{uj) = Q{4>ui) and q'{uj) = —q{uj)). (1.4) 

.2 doj 

Using these conservations laws, if (j)^ is a strict local minimizer up to constant phase e'® of E 
under the constraint | 1 w| 1 l 2 = ||</|cj||l 2 , then it has been shown that is orbitally stable, see 

[8, 31, 32, 80]. We are interested on standing waves not covered by this classical result. We will use 
the following notion. 

Definition 1.2. A bound state is not trapped by the energy if for any e > 0 there exists an 

Ue satisfying - Ue\\m < e, I|u£||l 2 = ||^<^||l 2 and E{ue) < E{4>^). 

The results [8, 31, 32, 80] do not cover the case when is a local but not strict minimizer. 
They also leave unsolved the case when (j);^ is not trapped by the energy, which we will discuss here. 

In order to study the stability, it is natural to consider the linearized operator C^i of (see 
(2.7) for the explicit form). Indeed, if C^i has unstable modes (spectrum with positive real part), 
then is orbitally unstable. (Even though this may look trivial, it is not, and it was proved 

rather recently first in 2D by [56] and later in general by [28]). Classical papers proving orbital 
instability of solitary waves by first proving their spectral instability are [30, 32, 37, 38]. 

In the case of ground states, except for the degenerate cases when q'{u)) = 0, only the above 
two cases (trapped by the energy or linearly unstable) occur. That is, if q'{uj) < 0 then has 
an unstable mode while if q'{(jj) > 0 it is trapped by the energy, [31, 32]. For the degenerate case 
q'{uj) = 0, see [9, 49, 61]. 

Excited state are usually not trapped by the energy and furthermore there are cases when 
Cui has no unstable modes. For example, if —A + V has n simple negative eigenvalues {ei < 
62 < ••• < en(< 0 )} then, if 26 ^ < Ci for some j > 2 , the excited states bifurcating form are 
not trapped by the energy and are spectrally stable. Even if spectrally stable they are orbitally 
unstable, [18]. Prior to [18] no systematic proof of this orbital instability was available. The series 
[26, 27, 34, 60, 71, 72, 73, 74, 75, 76], which stemmed from [7, 70], was able to treat only case 
262 > ei where, generically, excited states are spectrally unstable, see [20, 34, 60, 76]. 

Vortices (1.2) of (1.1) in the important pure power case /3 (|u|^)m = —\u\p~^u have been consid¬ 
ered in [57, 58, 59] which have various instability results, always by proving first spectral instability. 
Another important example is given by [63, 65, 77] for vortices for the cubic-quintic nonlinearity 

iut = —Alt — (|u|^ — |m|^)m, u(0, x) = uo{x), {t, x) € R x R^, (1.5) 

(for a review paper see also [10], for work on spinning solitons in 3D see [55]; see also [5]). [63, 77] 
show numerically for some values |?7i| > 1 the existence of a critical value Wcr such that for co < lo cr 
the vortices are spectrally unstable and for lo > LOcr are spectrally stable. In the simulations in [63] 
the spectrally stable vortices for lo > LOcr appear stable while in [77] for m = 3 appear to slowly 
develop instabilities. This latter observation appears consistent with the more recent numerical 
observations in [22, 43], in turn based on the instability theory in [14]. The theory in [14] is centered 
on the notion of Krein signature, which we introduce in Lemma 2.4 (although the standard definition 
is in the proof of Lemma 7.1). In this paper, we generalize [14] by using a simple idea from [18]. We 
will need the following notion. 

Definition 1.3 (Conditional asymptotic stability). We say that a bound state is conditionally 

asymptotically stable if there exist constants eg > 0 and Cq > 0 s.t. if u G C''([0, oo), H^) is a solution 
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of (1.1) with supoginf^gK ll“(^) ~ < e < eg then there exist w+ € O, 0 G and 

h+ G with ||/i+||//i + |a;+ — a;| < Cge such that 

lim \\u{t) - - e‘*^h+||//i = 0 . ( 1 . 6 ) 

t—>-oo 

Then we have the following orbital instability result. 

Theorem 1.4. Consider a bound state e^‘^*<f)uj and suppose that Hypotheses (H1)-(H5) in Section 2 
below are satisfied. Then if is not trapped by the energy and is conditionally asymptotically 

stable it is also orbitally unstable. 

Proof. Here the key hypothesis is (H5), i.e. q'{uj) 7 ^ 0 and hence either q'{uj) > 0 for all uj or 
q'{uj) < 0 for all ui. We prove the theorem by contradiction. 

Assume that the statement is false and that there is a wi G O s.t. the standing wave is 

orbitally stable. Since is not trapped by the energy there are initial data Mg arbitrarily close 

to such that E{uq) < and Q{uff) = Q{(j)uji). We can apply the Conditional asymptotic 

stability, Definition 1.3, and conclude by the conservation of the mass that 

q{uJi) = Q{u{t)) = lim Q{u{t)) = q{uj+) + 2 "i||h +|||2 > q{uj+). 

t—¥ao 

Similarly, by the conservation of the energy we have 

> Eiuo) = lim Eiuit)) = E{ct>^^) + 2-^\\Vh+\\l, > (1.7) 

t—foo 

By VE{(I>^) = —uj'VQ{4>uj) we have ^E{(j)^) = —uiq'{uj) (recall the notation q{uj) := Q{(j)u)) in (1-4)). 
On the other hand, by (H5) we know that q'{uj) has a fixed sign. So, since w > 0, both uj E{(j)i^) 
and UJ —>■ q{uj) are strictly monotonic, one increasing and the other decreasing. This means that 
<z(wi) > (z(w+) implies that E{(j)^f) < E{(j)^^^). But this contradicts E{(j)^^) > i?((/)^_^) in (1.7). This 
contradiction shows that is not orbitally stable for any wi G O. □ 

Remark 1.5. Notice that (after modifying the assumption (H4)) Theorem 1.4 holds for arbitrary 
dimension and we do not require that the nonlinear bound states are vortices. 

Remark 1.6. We do not claim that conditional asymptotic stability is necessary for the instability 
of excited states. 

Remark 1.7. For discrete NTS there are examples of bound states which are not trapped by the 
energy, are orbitally stable, but are not asymptotically stable, see [50, 51]. 

We now start the discussion on how to apply Theorem 1.4. We first discuss a sufficient condition 
for the standing waves to be not trapped by the energy. One natural way is to look at the Taylor 
expansion of the energy, which is often used in the study of orbital stability. Set Si^iu) = E{u) + 
ujQ{u) (notice that if Q{u) = Q{v), E{u) > E{v) is equivalent to Suj{u) > Sui{v)). Then, since 
is a bound state if and only if VS'aj(^cj) = 0 , we have S'tj(^!)tj+u) = Sui{4>lj) + ^ {y^Sui{(l)uj)v,v')+o{v‘^). 
Roughly speaking, we have one constraint Q((jju> + r") = QifjJui) which may eliminate at most one 
negative direction of V^S^{(jj^). Therefore, if has more than two negative eigenvalues, 

(j)^ is not trapped by the energy. Further, for the case has one negative eigenvalue, the 

trapping/nontrapping can be determined by the sign of q'{uj). That is, if q'{uj) > (<)0 then 
is trapped (not trapped) by the energy. 

Another viewpoint which we adopt in this paper is to study the linearized operator = 
given in (2.7), where J = f Here, we have extended V^Suii4>uj) to be a matrix 


3 


to make it self-adjont, see section 2. Since has the form J(—A+u;)+“rapidly decaying potential”, 
we see that \ae{Lui) = (—oo,—w] U [a;,oo). Here ae{Lui) is the essential spectrum of Lui- Since we 
are only interested in the case that is spectrally stable we assume C iK. 

In general no eigenvalues are expected in the interior of the essential spectrum. References [2, 54] 
have computational proofs about the absence of embedded eigenvalues when tpLj is a ground state 
for some equations. [63] does not discuss explicitly the issue of embedded eigenvalues but they seem 
to be excluded, see Remark 7.10 below for some further comments. [20] proves that in some generic 
sense embedded eigenvalues do not exist because they are unstable, see also [30, 76]. The proof is 
similar to the discussion of the well established instability of embedded eigenvalues in the case of 
self-adjoint operators, cfr. [40]. However it is not clear whether taking for example V generic in the 
cubic-quintic NTS would make generic in the sense of [20]. We think that embedded eigenvalues 
(with some negative Krein signature) can exist, but that they are unstable (we conjecture the non 
existence of embedded eigenvalues of positive Krein signature). Hence the assumption of absence of 
embedded eigenvalues seems reasonable. 

Similarly, the edge of the essential spectrum is shown in [19, 78] generically to be neither an eigenvalue 
nor resonance. Again, we did not check whether taking V generic makes the associated to the 
cubic-quintic NTS generic in that sense, although this is even more likely than for the issue of 
embedded eigenvalues. In [63] there is no discussion about the edge of the essential spectrum where 
the presence of a resonance or an eigenvalue would affect the computations discussed in pp. 371-372. 
We think that assuming absence of eigenvalue or resonance at the edge is reasonable. 

We will further assume that the kernel is spanned by and all nonzero eigenvalues (which are 
assumed to be finitely many but in fact this can be proved to be the case, see the comments 
under (Hll) in Sect. 2) have the same algebraic and geometric dimensions which are finite (in the 
presence of a nonzero imaginary eigenvalue whose two dimensions differ we would have a sort of linear 
instability, see Remark 7.10 for further comments). Now, let ^ be the eigenfunction of C^i associated 
to iA with A > 0 (for eigenfunctions of iA with A < 0, we have the symmetry. Thus it suffices 
to consider the positive case). Then, the “energy” of ^ is = iAH(^,^), 

where il = •, •) (see (2.6)). It is known that we can normalize ^ s.t. H(^, ^) = is for s € {1,-1}. 

Therefore, if s = 1, ^ have a negative energy and if s = — 1, ^ have a positive energy. Now, set 
n{\7‘^Suj{4>ui)) be the number of negative eigenvalues of V^Sui{4>ui)- It is known that n{V^Sui{4>u:)) 
can be represented as 


= p{q') + 2N- + A, + 2A„ 

where p{q') = 1 if g' > 0 and 0 if g' < 0, Ni, Nc are the number of real, imaginary and complex 
eigenvalues of (See [20, 79]). In our situation, Ni = Nc = 0. Therefore, if N~ = 0, (j)uj will be 
trapped by the energy and if N~ > 1, the is not trapped by the energy. Thus, by the above 
discussion, the trapping/nontrapping of the energy is determined by the existence of eigenvalue of 
with negative energy. Here, we will give a direct proof of the nontrapping when there exists a 
negative energy eigenvalue. 

Proposition 1.8. Assume hypothesis (Hl)-(Hll) and (HU;.) in Section 2. Then the standing waves 
are not trapped by the energy. 

Proof. Here (H14) is the key hypothesis which is roughly stating that for any w there is one eigenvalue 
iAj £ M+ with negative Krein signature (i.e. the corresponding eigenfunction satisfies H(^j , fj) = 
!)• 
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Fix uJi € O and choose a(e) to satisfy Q ((1 — a(e))i^i^^ + e(^j(wi) +fy(wi))) = Q{(j)u}^)- Since 
e Ng{Cl,J^ (see (2.9)), we see that 

Q ((1 - a(e))(()tji + e{^j{LOi) + ^j(wi))) = (1 - a{e)YQ {4 >ujt) + ((fy(wi) + fy(wi))) ■ 

Thus, we can conclude Qf(e) ^ e^. Consequently, we have 

E {{I - a{e))(l>^^ +e(^j(wi)+fy(wi))) - E{(l)^,) 

= 5'^i ((1 - a(e))(?i,^i +e(^j(a;i)+ fj (wi))) - 5^1 (<(>., 1 ) 

= y (fy(wi) + fy(wi)) ,fy(wi) + fy(wi)) +o(e2). 

Finally, since Cuji = we have 

(V2S'^^(,/>^J (^^(wi) +^j(a;i)) ,Cj(wi) +fy(wi)) = ^ (Cj(wi) + fy(wfy) , Cfywi) + fy(wi)) 

= iAjn(^j (cai), (lci)) iAjn(ify(cci), (cui)) = 2Aj < 0. 

Therefore, we have the conclusion. □ 

We now turn our attention to the conditional asymptotic stability. Its proof is very close to 
the proof of asymptotic stability of ground states in [12, 13] and we will need assumptions for the 
spectrum of Cuj discussed above (i.e. absence of embedded eigenvalues, ±iw are neither eigenvalues 
nor resonance, the generalized kernel is 2 dimensional, all nonzero eigenvalues have nonzero energy, 
see section 2). We will consider only the 2D case. 

Theorem 1.9. Assume hypothesis (H1)-(H13) in Section 2. Then, we have the conditional asymp¬ 
totic stability. 

The main example for the theory developed in this paper that we have in mind is the cubic 
quintic NTS (1.5) and specifically perturbations obtained adding a radial potential F(|a;|)M. As 
we mentioned above, in [63] it is shown that for any m = 1,2,3,4, 5 there is a critical value ujcr 
such that for uj < uJcr the vortices are spectrally unstable and for uj > uJcr are spectrally stable. 
[63] proves numerically that the linearization has an eigenvalue of negative Krein signature, i.e. 
our hypothesis (Hid) is satisfied. Indeed in [63] it is shown numerically that for (1.5) the spectral 
instability for w < Wcr occurs because a pair of eigenvalues on iR coalesce as w —>■ a;+ and then 
for Lo < ujcr form two eigenvalues which exit iR in opposite directions. This can happen only if the 
eigenvalues for uj > Ucr do not have the same Krein signature. Thus for lo close to and larger than 
LOcr at least one imaginary eigenvalue has negative signature. See below in Lemma 7.1 for a more 
precise discussion. Our results do not apply directly to the 2D cubic quintic NTS (1.5) because of 
its translation invariance. However, as we show in Section 7, when we add to (1.5) a small linear 
potential, then we obtain an equation which satisfies hypothesis (H14) and which is not translation 
invariant. When we take the small potential with a nondegenerate minimum at the origin, then for 
UJ > ujcr we obtain spectrally stable vortices. Our hypotheses (H1)-(H7) are either obvious or we 
know they are true as a consequence of the numerical experiments in [63]. The hypotheses (H9)~ 
(H13), while probably generically true, ought to be checked numerically. We will say more later 
about them, especially (H13), the most delicate and least analyzed. The conclusion that we can 
draw is that, assuming that indeed (H9)-(H13) are true, then for uj close to and larger than Ucr the 
vortices of appropriate perturbations of the cubic quintic NLS (1.5) are not trapped by the energy 
and are conditionally asymptotically stable. And hence by Theorem 1.4 they are orbitally unstable. 
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We have seen that Theorem 1.4 is a simple consequence of non trapping by energy (by Propo¬ 
sition 1.8 a consequence of the existence of one eigenvalue of negative Krein signature) and of the 
Conditional asymptotic stability. The latter follows from (H1)-(H13) by Theorem 1.9. A more 
precise formulation of Theorem 1.9 is in Proposition 2.5. The Conditional asymptotic stability, like 
the asymptotic stability of ground states proved in [13], is due to a mechanism of loss of energy of 
discrete modes related to the Fermi golden rule (FGR) and to linear scattering of the continuous 
modes. 

The FGR was first introduced in [7, 67, 70] and we will discuss it shortly. But first we discuss 
linear scattering, which is based on a number of results on the group 6*"^“ associated to the lin¬ 
earization of the NFS at the vortex. The results needed here are quoted (here we focus on 2D) from 
[21] and require that Lui should not have eigenvalues (and resonances) in the essential spectrum, see 
hypotheses (H9)-(H10) below. The results in [21] allow to say that, for all practical purposes, the 
restriction of on the continuous spectrum part, behaves like restricted to its contin¬ 

uous spectrum part. (H9)-(H10) are probably generically true, but nonetheless ought to be either 
proved or checked numerically on any given example. The fact that we cannot treat translation (the 
asymptotic stability result on moving solitons in [15] has not been proved in dimensions 1 and 2) 
depends on the specific way in which scattering of continuous modes is proved. Probably there is a 
simpler and more robust way to prove dispersion using virial inequalities and the theory of Martel 
and Merle [52, 53]. This should require fewer hypotheses and lead to similar results. This approach 
has been very successful in the context of KdV equations where it has improved the result by Pego 
and Weinstein [64]. The theory is proving successful also in different contexts, see [6, 11, 29, 47, 42]. 
However here we follow our standard approach and we use material in [21]. 

We now turn to the FGR. It can be easily seen under appropriate coordinate systems. It leads 
to nonlinear interactions between discrete an continuous modes of the NTS which are responsible 
for energy leaking out from the discrete modes. 

When we discuss the FGR we need to separate two distinct issues, as we will see with a simple 
example below. One issue is the fact that certain coefficients of the system have a 2nd power 
structure. This has been proved in [13]. See also [12, 18] for generalizations and references. A 
separate issue, is whether or not these 2nd powers, which are non negative, are also strictly positive. 
There might be cases when this is not true, but in general we expect that they are strictly positive. 
We do not have the expertise to run numerical tests, but a simple model might clarify this point 
(for other examples of FGR see also the survey [82]). 

For {z,h) € C X iJ^(R^,C), consider the Hamiltonian 

'H{z,h) = + \\Vh\\'j ^2 + \z\'^z ( G{x)h{x)dx + \z\^z f G{x)h{x)dx, (1.8) 

where G is a C-valued Schwartz function and the symplectic form idz A dz + 2 {idh, dh) where 

(/, 5 )=Re/' f{x)g{x)dx= ( {fi{x)gi{x) + f 2 {x)g 2 {x))dx iov f,g ^ C{=^^) . (1.9) 

Jw Jk? 

Then we have the Hamiltonian system 

ih = —Ah + \z\‘^zG, 

ii = z + 2|zp j h{x)G{x)dx + z^ f h{x)G{x)dx 

The solution of the linearized equation around (0,0) is {z,h) = {e~'*z{0),e'*^h{0)). Therefore, 
at the linear level, we do not see the asymptotic stability of (0,0). In the following, we sketch 


( 1 . 10 ) 

( 1 . 11 ) 
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heuristically why the equilibrium (0,0) is asymptotically stable thanks to the FGR and scattering 
of the continuous mode. First, if we assume z{t) = e“‘‘ 2 ;( 0 ), then — \z\^zR^^{l)G solves (1.10) (for 
the use of i?l^(l) = lim£_,.o+ i?_A(l +ie) see the remark on p. 30 [70]) . Therefore, it is reasonable 
to set 

h = -\z\'^zR+_^{l)G + g (1.12) 

and think 5 is a remainder. In fact (1.12) is a normal form transformation intended to eliminate 
the term \z\'^zG from the r.h.s. in (1.10). g satisfies an analogous equation as h, but with a higher 
degree polynomial in (z, z), and so it is smaller than h (this goes back to [7, 70] ) 

When we substitute (1.12) in (1.11) and we ignore g, we get 


iz = z-2\z\*z [ GRtA{^)Gdx-\z\^z [ GRZ^{l)Gdx. 

We recall that i?^^(A) = P.V.{—A — A)“^ ± i 7 rd(—A — A) for any A > 0, where S is the Dirac delta 
function and P.V.x~^ is the Cauchy principal value. These can be given sense using the Fourier 
transform. Multiplying by z and taking imaginary part we get 

^|zp = —27rc|z|® with c = f GS{—A — l)Gdx>0. (1-13) 

dt 7 r 2 

We conclude that c > 0 by the following formula, see ch.2 [23]: 

f \G{OfdaiO. (1.14) 


This is the 2nd power structure discussed above. In the context of the study of ground states of the 
NTS, the analogue of this formula has been proved in [13], see formulas (6.13)-(6.14) and later in 
this section. 

The next step is to ask whether in (1.14) we have not only c > 0 but rather the strict inequality 


c > 0. Obviously, for generic G we have G 


151=1 


^ 0, and so c > 0. Numerical computations for 


random choices of G would yield convincingly this fact. The strict inequality c > 0 and the equation 
of z in (1.13) yield the explicit formula 


W)? 


k(o)P 

(1 + 47rc|z(0)|^t)^ 


On the other hand, h will scatter by linear mechanisms, thanks also to the fact that the forcing 
term |zpzG in (1.10) is in for time. (Combining this fact with Kato smoothing estimates snch 
as Lemmas A.3 and A.4, we can show h has finite Strichartz norm, which implies the scattering). 
Proceeding differently and as a reference for later comments, one could integrate (1-13) and write 

|z(t)|2 + 2^c f |z(s)|®ds = |z(0)|2. (1.15) 

Jo 

Our FGR hypothesis, stated explicitly in (6.15), is the same as assuming c > 0 in model (1.8). The 
analogue of c > 0 instead is rigorously proved in (6.14). Numerical computations are likely to prove 
(6.15) true for generic equations exactly in the same way they would show that c > 0 in the above 
model. 
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We now discuss the FGR, still heuristically, for a model closer to the one necessary to examine 
equation (1.1). The discussion is more complicated than for model (1.8), but will yield similar 
conclusions. 

First of all, by a Noetherian reduction of coordinates related to the 17(1) invariance of the NFS, 
we will see that we reduce to an effective Hamiltonian of the form, for appropriate finite sums, 

n 

"^(2,/i) =-^ SjAjjzjP + ((-A + w + 9J)/i,cri/i) + ^ (1.16) 

j=l 


where: (z,h) € C" x i7^(]R^;C^) with h = aih; the inner product {f,g) = J *f(x)g(x)dx is a 
bilinear map; iAj G i]R+ are eigenvalues of the linearization A = (Ai,..., An); fU(a:) is smooth in 
X, rapidly convergent to 0 as x —oo and is for every x a self adjoint 2x2 matrix; w > 0 and in 
the application we are thinking it is in O of (H4); /x, z/ are multi-index such as g = (/ii, • • • , /x„) and 
z^ = Zi^ ■ ■ ■ z^" and similar for z''\ the second (finite) sum is taken for multi-indices /i, v] Gfiu{x) is 
Schwartz class in x with values in 2 components vectors; 


0-1 




(1.17) 


the number —Sj is the Krein signature of each iAj, with at least one Sj = 1 (for the connection 
between Krein signature and number Sj see in the proof of Lemma 7.1). 

Here h = aih, the Hamiltonian 'H{z,h) is real valued and so = —uiGufj.- We consider the 
symplectic form 

n 

i Sj dzj A tlzj -I- i{dh, a^aidh). 

The equations are then of the form 


ih = ICh + z“ 2 ^Ga,fl for/C = (T 3 (—A-|-a;-b 03) 

|A.(a-/ 3 )|>a; 


iSjZj =-SjXjZj + Vj— —(G^y,cr3(Ti/i). 


(1.18) 




Setting h = z°‘z^R^{X ■ (/3 — a))Gap + g like in (1.12), substituting and ignoring g (as we did 

earlier) we get 


\\-{ol-I 5 )\>u} 


As we will see, we can ignore the terms where X - — v) ^ X • [a — P). Furthermore, up to smaller 

terms that we ignore, we have 

iSjZj =-SjXjZjUj^—{Gou,(Js(JiRK{-L)GM). 

L>uj X v—X-a—L ^ 

Let us write formally R'^{—L) = P-V.-j^^ -I- 'nrd{IC -I- L) (there is a distorted Fourier transform that 






allows to make sense of this). Then, using Gqi/ = —criG^o, 

=-TT^ LRe/ ^ z^Go.,c73C7iS(/C + L) ^ z“G„o\ 

j=l L>oj \x-l'—L Xa—L / 

= -TT L{Gl{z), g3(5(/C + L)Gl{,z)) where Gl{,z) := ^“Gao- 

L>uj X-Oi—L 

Furthermore, there exists G^^\z) s.t. 


(1.19) 


{GL[z),a35{K. + L)GL{z)) = {Gf{z), ^(-A + w + L) 0 Gf{z)). 

\ 0 -(5(A-a; + L)/ 

Observe now that (5(A — uj + L) = (5(—A — {L — w)) and that integrating in (1.19) we get 
-TT^] T f {{G^°\z{s)))^,S{-A - {L - uj)){G^°\z{s)))^))ds 

j=l L>lo 

n 

i=i 


( 1 . 20 ) 


Notice that inside the integral we have a similar 2nd power structure (for each L we have a positive 
quadratic form in the vector {z'^)\.a=L ) to the (1-14) we found in model (1.8). In particular, if n = 1 
in (1.16), then the time integral in (1.20) will be similar (possibly with different power) to the time 
integral in (1.15). So formula (1.20) is a generalization of (1.15). When the Krein signatures are 
all positive, that is if = — 1, then (1.20) can be used to prove z{t) 0 (and so the asymptotic 
stability of the standing wave). Indeed, starting from a |2:(0)| small, (1.20) is telling us that also 
\z{t)\ remains small. Furthermore, the fact that the integrals remain bounded as t —>• oo can be used, 
along with the fact that \z{t)\ remains bounded (which can be seen by the 2nd equation in (1.18)), 
to show that z(t) q. However, to make this rigorous we need to have something analogous to 
the inequality c > 0, which is exactly the meaning of (6.15). 

If there are negative Krein signatures, and so some Sj = I, obviously the proof of z(t) q 
breaks down since in the l.h.s. of (1-20) there are terms with different signs whose size could become 
large even if the sum is small. 

However, if we know that a solution u{t) remains close to a standing wave, and consequently 
the corresponding |z(t)| remains small, then (1.20) allows to prove z{t) q because we know that 
the sum '^j SjXj\zj(t)\‘^ remains small. This allows to control the integrals and, with hypothesis 

(H13), that is (6.15), to conclude z{t) 0. 

The validity of (6.15) (or of more explicit formulations of the formula) is an open question. 
Proofs related to special situations are in [3, 46, 1]. It is fair to expect that numerical experiments 
will confirm strict positivity for most examples, like for the c in (1.14). Some numerical verifications 
are in [22, 43], but there is room for more systematic studies. These are absent in the literature not 
because some intrinsic difficulty, but simply because the ideas in [13, 14] are not well known. 

Having given a general overview of the main concepts and results of this paper, we list now 
the content of the remaining section of the paper. In Section 2 we state hypotheses and give in 
Proposition 2.5 a statement that is more detailed than Theorem 1.9. In Sections 3-6 we describe 
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the proof of Proposition 2.5. The proof is basically the same of the proofs of asymptotic stability of 
ground states of the NTS in [12, 13]. In Section 3, after introducing a natural system of coordinates 
related to the modulation we state in Proposition 3.4 a result on Darboux coordinates, whose proof 
is in [12] and which is a key step to the subsequent search of an effective Hamiltonian. We will only 
give a sketch here. The expansion in these new coordinates of the functional K{u), defined in (3.14), 
is again given without proof since the proof is in [12]. In Section 4 we complexify L^(R^,]R^) and 
consider the spectral decomposition of the linearization C^i in L^(R.^,C^). This produces discrete 
coordinates z = G C" and a continuous coordinate / G We then state the 

expansion of the functional IC(u) in these coordinates, with the proofs in [12]. Next step is the search 
of an effective Hamiltonian by means of a Birkhoff normal forms argument. This is accomplished in in 
Proposition 4.3, which we give only a sketch of the proof, for which we refer again to [12]. Proposition 
5.1 contains Strichartz and smoothing estimates satisfied by /, and the statement of z(t) Q. We 
then show that Proposition 5.1 implies Proposition 2.5. By an elementary continuation argument 
Proposition 5.1 is in turn a consequence of Proposition 5.2, whose proof is in Section 6. In consists 
first in Strichartz and smoothing estimates for /. We then split / like in (1.12) as a term dependent 
only of z, which is the part of / which really affects the z’s, and a g which is smaller, does not affect 
z substantially and can be treated as a reminder term. The estimates on / and on g are the same 
of [21]. Finally in Section 6 we return to the Fermi golden rule, explaining why z{t) 0. All the 
estimates are proved in the literature, for example in [16]. Therefore, we limit ourselves at describing 
the structure of the argument. However we give a sketch of the proof for some important theorem 
(especially Darboux theorem and Birkhoff normal forms arguments) for reader’s convenience. 

In Section 7 we discuss the cubic quintic equation (1.5). We discuss how starting from the 
numerical observations in [63] it satisfies hypothesis (H14) for values w > LOcr close to ujcr- Since our 
theory does not apply to translation invariant equations like (1.5) we show that by adding a small 
radial potential with a non degenerate local minimum at 0 produces spectrally stable vortices which 
still satisfy (H14) because their linearization is a small perturbation of that of (1.5). At the end of 
Section 7 we also discuss the status of the other of the other hypotheses for equation (1.5) perturbed 
by adding a liner potential. Some of them follow from the computations in [63], the others ought to 
be checked numerically and in out opinion are plausibly true. 


2 Hypotheses and statements 

To begin with, assume the following hypotheses. 


(HI) 

(H2) 


;9(0) = 0, /3 G C'°°(K,R). 

There exists a p G M such that for every k > 0 there is 


dv^ 


Piv^) 


< Ck\v 


p—k — 1 


a fixed Ck with 
if |u| > 1. 


( 2 . 1 ) 


(H3) V(x) is smooth, non zero, real valued, and for any multi index a there are Ca > 0 and Oa > 0 
such that |i9“H(x)| < 

For n > 1 and AT = M, C then = S„(]R^, AT^) is the Banach space with 

IkllL := + ||5“u|li2(R2)) < oo. ^2.2) 

|Q;|<n 
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We set So = We define S* by ||m||ej := ||(1 — A + |xp) 2 u||i 2 < oo for t G R. For t G N 

the two definitions are equivalent, see [15]. 

(H4) There exists an open interval O C K+ such that 


Au — Vu — uju + P{\u\'^)u = 0 for a; G (2.3) 

admits a function ui ^ (p^, which for any k is in Sfc(]R^, C)). We also assume that 

'Pui{x) = like in ( 1 . 2 ). 

(H5) We have for q{u}) := 2 “ 1 ||(^,^||| 2 (r 2 ), see (1.4), 

g'(w) ^ 0 for all to G O. (2.4) 


Remark 2.1. Notice that by a standard bootstrapping argument we can relax hypothesis (H4) by 
only asking that uj pui be in (7^(0, iJ^(]R^, C)). 

We identify C = setting wi = Re w and W 2 = Im w for w € C . In particular we identify the 
imaginary unit i with — J where 


J = 




(2.5) 


and the bound state (pui with 


Re0(, 
Im A 


We consider the strong symplectic form defined 


n{X,Y) = {J-^X,Y). 


( 2 . 6 ) 


Definition 2.2. We denote by ( , ) also the bilinear form in L^(R^,C^) obtained extending (1.9). 
We extend Q to L^(R^,C^) as a bilinear form. 

For F G C^CUjR) with U an open subset of i7^(R^,R^), the gradient VF(it) is defined by 
(VF(m),X) = dF{u)X, with dF{u) the Frechet derivative at u. If F G C^(U,R) it remains de¬ 
fined V'^F{u) G (^^(U, i?(i7^(R^, R^), iJ~^(R^, R^))) (with il(X, Y) the space of R-linear bounded 
operators from a Banach space X to another a Banach space Y). 

The (puj are constrained critical points of E with associated Lagrange multiplier —uj so that 
XE{(pi^) = —u!XQ{(pi^). The linearization of the NTS at (puj is 

:= J{V^E{(P^) + 00 ) = J(-A + u: + V)FJVu. with (2.7) 

V P'{\(pu,\‘^){b{.epu,){lm.(p^) /3(|(/ii^p)-I-2/3'(|0^P)(Im0^)2y ■ 


Since there is a natural identification L^(R^, C^) = L^(R^, R^)®rC and i7^(R^, C^) = i7^(R^, R^)( 8 )r 
C, the operator extends naturally in a C linear operator in L^(R^,C^) with domain 77^(R^,C^) 
simply starting from formulas £(^(n (8 )r z) = {L^^v) 0r z. 

Starting from v ( 8 )r z := n( 8 )RZ a complex conjugation can be defined in F^(R^,R^) ( 8 )rC, which 
should not be confused with the complex conjugation in the initial L^(R^,C). Using this complex 
conjugation we consider in L^(R^,C) the hermitian form (/, 5 ), where we recall that ( , ). is a 
bilinear form. 

Notice that if Cu)^ = z^ with z G C, then applying this complex conjugation we obtain = z^. 
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Notice also that as operators in L^(K^,]R^) we have JC^j = —J and C^jJ = —. This extends 
also for the operators in L^(]R^,C^). From this we conclude that (j{Cui) is symmetric with respect 
to the coordinate axes. 

We consider only standing waves which are spectrally stable. Specifically, we will assume 
(H6) a{£ui) C iK for all uj G O. 

Since (H4) implies that is exponentially decreasing to 0 as a: —^ oo, we know that always 

for the essential spectrum we have i(Je{Cu)) = (—oo, —w] U [aj,oo). Hence (H6) is all about the set 
of eigenvalues ap{Cuj). We want our to be also linearly stable. This is somewhat ambiguous 
as in principle it should mean that ||e*'^“ ||i 2^^2 is bounded, which is never true since is not 
skew-adjoint and has a nontrivial Jordan block at 0. By linear stability we mean (H6) and the 
following two additional hypotheses. 

(H7) For Ng{L) := ker(L'’) we have 

ker£^ = Span{ and Ng{C^) = Span{ (2.8) 

(H8) For any eigenvalue e £ crp(£ij)\{0} the algebraic and geometric dimensions coincide. 

Notice that Ng{C^) for (1.5) has been computed in [63]. For what happens when a potential is 
added to (1.5) breaking the translation invariance see Lemma 7.4. 

We assume the following hypotheses. 

(H9) There are no eigenvalues of C^i contained in Ue{Cuj)- 

(HIO) The points iiw are neither eigenvalues not resonances of i.e. if C^F = iiwF in a 
distributional sense for an F G L°°(R^), then F = 0. 

As we mentioned in the introductions oth conditions appear to be generically true. 

We have the symmetry J£^ = —£() J. Thus, (2.8) implies 

= Span{((),^, (2.9) 

We have the following beginning of Jordan blocks decomposition, where we use the hermitian from 
{f,g) to define the orthogonal spaces, 

L'^{m.^,C'^) = Ng{£^)®N^{C*J . (2.10) 

[21], as consequence of [41], proves that lAri(£.) 1 |l 2^^2 is bounded under (H6)-(H10). Appro¬ 
priate dispersive and Strichatz estimates can be proved for the restriction for the space 

L^(a;), see Lemma 2.3 below. 

We assume existence of non zero eigenvalues. 

(Hll) For any oj G O there is a number n > 1 and positive numbers 0 < Ai < A 2 < ... < An such that 
(7p{£uj) consists exactly of the numbers ±ej and 0, where we set 6 ^( 0 ;) := iAj(w). We assume 
that there are fixed integers no = 0 < rii < ... < = n such that \j = Xi exactly for i and j 

both in (n;,n;+i] for some I < Iq. In this case dimker(£tj — iAj(a;)) = n;+i — n;. We denote 

by Nj G N the number such that Aj- + 1 = inf{n G A/”: niXj G cre(£w)}. We set N = sup^ Nj. 
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Notice that in (Hll) we do not ask any more uniformity with respect to w, as in [13, 15, 16, 17]. 
We remark that the fact that the sum of all the algebraic dimensions of the eigenvalues of is 
finite can be proved from (HIO) and from the fact that each 4>u>{x) is in fact not only a Schwartz 
function in iS(M^,C) but converges exponentially to 0 as a; —>• oo, see [35]. The proof is standard, 
is similar to an argument in p.305 [66] and involves extending the resolvent beyond the resolvent 
set as a meromorphic function. Since we are in dimension 2, the discussion of what happens near 
iiw is more complicated than the 3D argument near 0 in [66], but nonetheless an accumulation of 
eigenvalues near iiw can be excluded using iiw. 

The following is a rather standard non-degeneracy hypothesis in the context of normal forms argu¬ 
ments. 

(H12) For distinct < ... < and fj, G with \fj,\ < 2N -t- 3, then 

/iiAji = 0 /X = 0 . 

It is plausible that (H12) is generically true. 

Next we assume the Fermi golden rule which we will state explicitly later and on which we commented 
at length in the Introduction. 

(H13) The Fermi golden rule Hypothesis (H13) in Section 6, see (6.15), holds. 

So far the hypotheses (H1)”(H13) are similar to the analogous ones in [13]. In [13] though 
the main result is that the standing waves are (asymptotically) stable, while here we want to prove 
instability, that is the opposite. So we need an hypothesis which will generate orbital instability. To 
obtain this hypothesis we consider the signature, or Krein signature, see [44, 45]. 

Recall the extensions of ( , ) and H in L^(R^,C^) made in Def. 2.2. Recall that ap{£uj) = 
ap{C^). By general argument we have the following result. 

Lemma 2.3. The following spectral decomposition remains determined: 

L^(R^,C^) = L^(£i^) 0 L^(£^) where LI{u}) = and for L = £^, £i^* (2-11) 

Ll{L) := Ng{L) 0 Ll{L) with Z^(£) := 0eGap(£,.)\{o} ker(£ - e). 

We denote by Pc{u}) the projection on L‘^{uj) associated to (2.11). 

The form D remains symplectic also in £^(£;^). The proof of the following lemma is elementary, 
see for example Lemma 5.2 [12]. 

Lemma 2.4. For any uj G O and corresponding n in (Hll) there are functions ^j(uj) G Sk for any 
k and j = 1,..., n such that the following facts hold. 

(1) G ker(£^ — iAj(a;)) for all j. 

(2) H(^j (w), ^fc(w)) = 0 for all j and k and = isjSjk with Sj G {1, —1}. 

□ 

In the case of ground states, that is when = tpui{\x\) with ifaiHx]) > 0 and to = 0, then 

Sj = — 1, see Lemma 2.7 [17]. Here, where to 0, we assume instead what follows. 

(H14) There exists at least one j s.t. Sj = 1. 

We have already discussed, and we will say more in Section 7, that it has been shown numerically 
that this hypothesis occurs for spectrally stable vortices of the cubic-quintic NLS (1.5). 

Theorem 1.9 is a consequence of the following proposition, which is a consequence of [12, 13, 
16, 21]. 
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Proposition 2.5. Let uji £ O and assume (H1)-(H13). Then there exist constants eo > 0 and Co 
s.t. if u G C°([0, cxa), is a solution of (1.1) with supt>Qinf^gR l|w(t) — e'"^ < e < co then 
there exist uj+ & O, 9 £ and h+ G with ||/i+||jji + |aj+ — wi] < Coe such that 

lim llM(t) - - e"*‘^h+\\m = 0. (2-12) 

t^CO 

It is possible to write u{t,x) = +A{t,x) + u{t,x) with |A(t,a:)| < CN{t){x)~-^ for any N, 

with limt_>oo C]\[{t) = 0, with limt_).oo w(t) = a;+, and such that for any admissible pair {q,p), i.e. 

l/q+l/p=l/2,q>2, (2.13) 

we have 

II^IIl«([o ,oo),W^'P) — Coe. (2.14) 

From Section 3 to Section 6 we focus on Proposition 2.5. There are various steps. We aim at 
showing that there exists an effective Hamiltonian of the form (1.16). This has to be found through 
a Birkhoff normal forms argument, see Theorem 11 [39]. In order to initiate the process we need to 
to find and appropriate system of Darboux coordinates. 


3 Modulation and Darboux coordinates 


Asymptotic (or conditional asymptotic) stability arguments require traditionally, since [69], the 
choice of appropriate modulation coordinates. Indeed, since we are discussing the stability of vortices, 
it is natural to express a solution u(t) which is close to a vortex as a sum of a vortex plus an error and 
to frame stability in terms of what happens to the error. This and more is what modulation aims to 
do. The hrst step to define precisely this vaguely stated aim is the following standard Modulation 
Lemma. 

Lemma 3.1 (Modulation Lemma). Fix n G Z, wi G C and where O is given 

in (Hf) and J is given in (2.5). Then there exists a neighborhood tin of 4*1 in S_ji(IR^,M^) and 
functions uj G C'^{Un,0) and H G C°°(Z4^, R) s.t. u;(4'i) = wi and 'd(4'i) = and s.t. Vu G Un 

u = + R) and R G A^^(£* ). (3.1) 

We give a sketch of the proof. See also Lemma 2.4 [12] or Lemma 2.2 [16] for a detailed proof. 

Proof. It suffices to apply the implicit function Theorem to 


Fi{'d,u},u) 



will be invertible because of (H5). □ 

The above Modulation Lemma is the starting point to find appropriate coordinates in the 
neighborhood of 4'i in iJ^(R^,C). Solutions u{t) starting close to d/i will admit a time dependent 
decomposition (3.1). If u{t) stays close to the orbit of d'l for all time and scatters to a vortex, this 
will be equivalent at showing that R(t) scatters to 0 as t —>■ oo and u}{t) for some a;+ G O. 

This will be proved working on the parameters in the r.h.s. of (3.1). 
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Equation (1.1) can be expressed as u = J\7E{u). The following discussion is standard, and is only 
sketched in order to give an overview of the use of the parameters By substituting the 

r.h.s. of (3.1) and using we obtain after standard computations 


(w - +R = C^R + N{R) - (w - ^)JR, (3.2) 

where N{R) contains terms which are quadratic or higher order for R small. Denote by PNg(oj) = 
PNg{c^) the projection on Ng[C^) related to (2.10). Because = q'{oj), we have 

PNg{oj)X = +n{Jcj)^,X)dM, yXGS'iR^C^). (3.3) 

Applying the projection 

P{u;)-.= l-PNg{uj). (3.4) 

to (3.2) we obtain for R the following equation 

R = £^R + P{u})N{R) - (w - ^)P{uj)JR. (3.5) 


Since, as is well known, the term (oj — i?) is higher order in R, we can think of (3.5) as a perturbation 
of R = CujR- It is natural now to look at the rest of the spectrum of C^i- The main difficulty is to 
show that the discrete components of R associated to the point spectrum of Cgj, which at a linear 
level want to oscillate like the e“**2(0) component of the linearization of (l.lO)-(l.ll), will lose 
their energy because of some friction originating from the nonlinear interaction with the continuous 
components of R. This effect will be captured by an argument similar to the Fermi golden rule 
discussed in the Introduction. For that argument to work we need to find an appropriate system of 
coordinates. 

Lemma 3.1 does not provide coordinates. We co back to the projection PNg{uj). We have 
n{PNgiu})X,Y) = Q{X,PNgiuj)Y). By (H4)-(H5) and (2.8) for S{R'^,K^) = nfe>oSfe(R^ A'^) the 
space of Schwartz functions and for 5'(R^, AT^) = Ufc<oFfc(K^, AT^) the space of tempered distribu¬ 
tions, we have 

PNgiio) G C°“(C>,B(5'(K^A:2),5(R^iG2))) for Ar = R,C. (3.6) 

For P{uj) defined as in (3.4) we have uj —>■ P{uj)P{uji) G C°°{0, B{I]k,Ek)) for any k. By (H4) we 
have P{uj) T’(wi) in the operator topology of B{Tjk,Yk)- Thus, writing 

P(c^)P(c^l) = (1 + (P(c^) - P(wi)))P(c^i), 

we see that there exists an a > 0 such that if jw — wil < a the map P{uj)P{uji) restricts into an 
isomorphism from fl Xk to A^(£*) fl Xk for any k > —n with Xk equal either to El^ or 

to Efe. Hence for k > —n the map 

R X {a; : la; - wil < a} X n X,) ^ A^, 

(d, a;, r) -)■ M = + P{uj)r), 

is for ||r||xfc < a a local C°° diffeomorphism in the image. Therefore, (d^uijr) in (3.7) provides an 
initial system of independent coordinates. 

If we consider the function Q = Q{u), the map {d,oj,r) — {'d,Q,r) is a local diffeomorphism 
because of the assumption (H5). Indeed, applying implicit function theorem to 

P 2 {Q,p,r,uj) = - PNg{uj)r) + p+{(l)ui- Pxg(w)r-,r) - Q, (3.8) 


15 


there exists uj{Q,p,r) which is a smooth function defined in the neighborhood of Oj 0) G 

K X R X S-n- Notice that F 2 (Q, Q(r), r, uj) = Q{cj)^ +P{oj)r) — Q. We have put an auxiliary variable 
p because if we directly put ui{Q,r) to be the implicit function of F 2 = + P{uj)r) — Q, then 

we will only able to define w in the neighborhood of 0) in R. x L^. Differentiating, (3.8) by 

Q, p, r, we have 

Dquj = -dpw = (3.9) 

(V.u;, X) = A-i ((P(cc)r, Pw, (cc)X) + {Pm^ (u;)r, X)) , (3.10) 

where A= F P{uj)r, d^cj)^ + d^P{uj)r). 

We now expand D by using the coordinates (z?, Q, r). Notice that 

X = duX = d^uXg + OquXq + (VrW, Xr) , 

where X^ = ddX^Xq = dQX and Xr = drX. Then, after some cancelations, we obtain 

n = — d-d A dQ + n{P{uj)dr, P{u;)dr) (3-11) 

+ A~^dQ A n{dujP{uj)r, P{ui)dr) + (VrW + dqUJr, dr) A D,{di^P{uj)r, P{uj)dr). 

Notice that the coordinates {d,Q,r) are not a system of Darboux coordinates for the symplectic 
form n. 

We now prepare some notations. 

• Let F be a Frechet differentiable function. Then, its hamiltonian vector field is defined by 
il{Xp, Y) = dF{Y) for any given vector Y. In particular, we have Xp = JVF. 

• For F, G two scalar valued functions, we set the Poisson bracket by {F, G} := dF{Xc)- 

• It Q has values in a given Banach space E and G is a scalar valued function, then we set 

{g,G} -dGiXc). 

In the coordinate system {d, Q, r) our NTS can be expressed as 

Q = {Q,E} = iXE)Q, d = {d,E} = iXE)^, r = {r,E} = iXE)r. (3.12) 

Further, comparing the coefficients of Yg in il{XE,Y) = dEY by (3.11), we have {Xe)q = 
dQXE = 0. Therefore, we have Q = 0. Notice that this shows that with coordinates {d,Q,r) we 
have achieved a reduction of order in the system, see [62] p. 412, effectively reducing to the variable 
r only. 

In the sequel we choose wq such that if uq is the initial value in (1.1), then 

Q{cf^,)=Q{uo). (3.13) 

We consider (recall < 7 ( 0 ;) = Q{4’uj)) 

K{u) : = E{u) - E ((()^o) + lo{u){Q{u) - ^(wo)). (3-14) 

Then, since (Xq)^ = —1, {Xq)q = 0 and (Xq)^. = 0, we see that (3.12) is equivalent to 

(5 = 0, d = {-d, Xj+w, r = {r,K}. 

See [12] Lemma 2.6 and Section 2.3, and it is important that Q{uq) = q{uJo)- 

In the sequel the changes of coordinates will differ from the identity transformation by pertur¬ 
bations that can be written in terms of the two classes of symbols which we introduce now. 
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Definition 3.2. For I an interval with 0 in the interior, A C RxRx j^)nS_„) a neighborhood 

of (g(u;i), 0,0), we say that S' G C""(/ x yf, R) is if there exists a C > 0 and a smaller 

neighborhood A' of of (g(wi),0, 0) s.t. 

\S{t,Q,g,r)\ < C'||r||^_^(||r||E_„ + |£»| + |Q - g(wi)|)‘ in / X A'. (3.15) 

We will write also S = or S = ^?! ^)- Given a function F : I x —>■ M for 

a neighborhood of in we say that F = if there is a function s.t. 

F(t, u) = Q, Q{r), r). 

We say S = ^n;'oo if i? = for all to. We say F = if we can take n arbitrarily large. 

If F = for any to, we set F = F^’F 

Definition 3.3. A T G C™'{I x A, Sn(R^, C^)), with I x A like above, is and we write as 
above T = or 1 = Q, g, r), if there exists a C > 0 and a smaller neighborhood A' of 

(Po,Po,0) s.t. 

||T(t,Q, £<,r)||E„ < G||r||^ 2 _„(ll’'lls_„ + Ifiil + IQ - <7(wi)|)* in / x A'. (3.16) 

We use notation T = T = and 1 = S®’f as above. As above, given a function T : I x 

Uuj^ —>■ S-n we write F = if there is a function s.t. T(t,u) = Sl;^^(t,Q,Q(r),r). 

Next we consider the following symplectic form: 

fig := —di? A dQ + n(dr , dr ). (3.17) 

This is how our symplectic fl form should look in appropriate coordinates, indeed in Section 3 [12] 
the following Darboux Theorem is proved. 

Proposition 3.4 (Darboux Theorem). There is a local dijjeomorphism F around (puii to L^(]R^,]R^) 
such that F*Q = Og and which in the (d^Q^r) coordinates is of the form 

d' = d + F°’^Q,Qir),r), Q'= Q , 

□ 

Remark 3.5. Notice that the idea of taking as fixed point rather than as in [12], is taken 
from [4]. The proof of Prop. 3.4 is unaffected. 

For the convenience of the readers we give a sketch of the proof here. 

Sketch of the proof of Proposition 3.4- To make a change of coordinate to convert the symplectic 
form O into fig we need three steps. First, we find a 1-form F s.t. 17 — Og = dF. Next, we solve 
ix‘^B = —r, where ixaj{y) ■= uj{X,Y) and 17s := — 17g). Finally, let J^s be the flow of 

■£ys = A’®(>’s)- Then, we have 

4- = y: [Cx^^s - d.ris) = y:idtx^n, - dr) = o. 

ds 

Thus, y := y^ gives us the desired transformation. This is a standard proof of Darboux theorem 
(see [39]). 

In our situation, we have to care about the regularity of the transformation (or in other words, error 
from the identity). Therefore, we need to compute F rather explicitly. 
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First, we seek F. It suffices to find some F satisfying F = B — Bq + dC, where 

2Bo = Qd-ff + n{r, dr), 

2B = niu,-) 

= Qdd + il{P{uj)r, dr) + P{u})dr) + + P(oj)r, dc^cjiui + duiP{ui)r)duj. 


It is elementary that dB = Q, dBo = Ho- We have 


2{B - Bo) =d {il{<j)^,P{u})r)) + ^}{P{uj)r, d^Pr)dQU]dQ 

+ r2(—PjVg {uj)r + Vl{P{uj)r, dujP{io)r)dpOjJr + n{P{uj)r, 9^P(u;)r) Vr-w, dr). 

Therefore, we can choose F = n(rj., dr) + TqdQ as 


2Tr = —PNg {i^)r + ^{P{(jj)r, duiP{u})r)dpUjJr + f2(P(a;)r, di^P{uj)r)y 
2rQ = n{P{uj)r,du,Pr)dQU}. 

Since PNgiuj)r = Ppj^{uj) (P(wo) — P{uj))r and jw — wqI ~ Q{,r), we see F^. = S^d _|_ 7 ^ 0-2 j^, 
rQ = pod. 

Next, we solve ix‘^s = —T. Since s(II(T’'*,-) — f2o(T’'’,-)) can be handled as a perturbation, the 
main part of the equation will be r2o(-T'*,-) = —F. Therefore, we have Xq = 0, Xg = 

X^ = _j_ 7 ^ 0,2 Finally, solving = A’®(3^s), we have the conclusion. □ 

It is well known that normal forms processes are based on Taylor expansions of the Hamiltonian, 
see [39]. So we need an expansion of the functional K defined in (3.14) in terms of the Darboux 
coordinates. This is provided by the following, proved in Lemma 4.3 [12]. 

Lemma 3.6. Consider an integer L G Af s.t. L > p for p the exponent (2.1) in hypothesis (H2). 
For any preassigned {k,m) and in the coordinates ((5,d,r) of (3.18), K admits the expansion 


d(uj) - d(uJo) - (uj - ujo)q{uio) + -il{CujP{ui)r, P{ui)r) + Q{r),r) + Ep{P{uj)r) + R", 

L-l „ 

R" :='^{Bd{Q,Q{r),r),{P{uj)r)‘^) + / BL{x,r{x),Q,Q{r),r){P{uj)r)^ {x)dx with: (3.19) 

d=2 


• d(w) = E{cj),^) + ujq[uj); 

• P2(<?(wo),0, 0) = 0; 

• {P{uj)r)‘^(x) represent d—products of components of P{uj)r; 

• Bd{-, Q, Q, r) G C'™(W_fc, Efc(R^, i?((R^)®‘^,R))) for 2 < d < 4 with U-k C a neighborhood 
of ('?(wi),0, 0) m R X R X (iV^(£* J n S_„); 

• for C G R^ and (Q, g, r) G U-k we have for i + j < m 

£*) ?’)|Ib(e® 3 _2j^(K2_B((K2)»t^R)) < Ci- (3.20) 


Remark 3.7. We have d{uj) — d{uJo) — {uj — u:o)q(uJo) = 0(uj — wq )^ = Tif'^{Q{r)) +77.^’^ (Qj Q{'>^)x)- 
In Lemma 4.3 [12] inequality (3.20) is stated for |C| < e for some small e > 0, but in fact in the 
proof is unnecessary, thanks to (H2). Notice also that L = 5 in [12], but a similar proof holds for 
our choice of L. 
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Sketch of the proof. Notice that here we are just expanding K{u) = Sa;{u) — Sa;{4>uio) such as 


K{4>ui + P{u})r) = P{u})r) + o (^{P{uj)rf^ - d{uJo) + (w - wo)g(wo), 


where we have used ySu,{(fuj) =0. □ 

If we expand P{uj)r = r + {P{uj) — P{uji))r we obtain what follows, see Lemma 4.4 [12]. 
Lemma 3.8. The expansion of K in Lemma 3.6 can be rewritten as follows,with similar notation: 


^ = '^k miQ^ <5W) + 2 ^Ll{c^,r, r) + Q{r),r) + Ep{r) + R', 

L-l 

R':='^{Bd{Q,Q{r),r),r‘^) + BL{x,r{x),Q,Q{r),r)r^ {x)dx . (3.21) 


d=2 


4 Spectral coordinates associated to 

Recall that r G A^^(£* JnL^(]R^,]R^). We consider the spectral decomposition of r in terms of Cui^: 

n n 

r = ^ Zj£,j{uJi) + + f where / G Lliuji) and also / G (4.1) 

i=i i=i 

This yields new coordinates r —>■ (z, /) to replace r. 

Correspondingly we have the expansion 

n 

Ll{dr ,dr ) = ^ i Sj dzj A dzj + fl{df , df ). (4.2) 


Equation r = {r, K} splits into 

iz,=s~K, f = {f,K}, (4.3) 

where we recall that Sj G {1,-1} and Sj = 1 for at least one j. 

We have reduced our NLS to system (4.3). Obviously, having replaced r with (z, /), we need to 
rewrite the expansion of K in Lemma 3.8 in terms of (z, /). This is done in Lemma 5.4 [12], which 
we quote. 

Lemma 4.1. In the coordinate system (z, /) near (0, 0) for any preassigned pair (k, m) we have an 
expansion 


4 

^ = '^k’!L(Q’ Qif)) + H 2 + ^ Rj + Qif), f) with what follows. 

l=-i 


(1) For Q = Q{f) 


iL' =-^sA|^,r+ ^ 7^l;°„(g,p)z^z^ + 2-l^l(£^J,/). 

J —1 2 


(4.4) 


(4.5) 
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(2) We have for g = Q{f): 


R- 1 = E E ^^r{JSl’°^{Q,g),f)- 

\ti+v\=2 \ii+v\ = l 

for N as in (H12), g = Qif), g^j.i,{Q,g) = Tl^kln^Q^s), Gf^i,{Q,g) = S°’^((5, p) and with the 
symmetries and = —G^y, we have 

2N4-1 2N 

Ro= Y. ^i = i E 

|/i+j/|=3 |/i+i^|=2 

^2 = with B 2 (( 7 (a;i), 0 ) = 0, 

where f‘^{x) represents schematically d—products of components of f; 


Rz 

Ra 


E E z^r{Jsl’YQ,z,eJ)jy, 


=2iV+2 
L-1 


=2iV+l 


V {BdiQ, z, g, /), /") + f Bl{x, fix),Q, z, Q{f)J)f^{x)dx + Ep{f), 

d=2 


where the B’s are like in Lemma 3.6. 


□ 

Now we start discussing about the normal forms argument. It will consist in eliminating as 
many terms as possible from Rj with j = —1, 0,1. 

We set, for n the number associated to wi in (Hll), 

e = (ei,...,en). (4.6) 

Some of the monomials in Rj with j = 0,1 cannot be eliminated because they are resonant, that is 
of the following type. 

Definition 4.2 (Normal Forms). A function Z{z, g, f) is in normal form ii Z = Zq + Zi where Zq 
and Z\ are hnite sums of the following type: 

Zi=i ^ z'^r(JG^,(p),/), (4.7) 

e-(p-v)^lTe{CuiB 


where for fixed fc,m G N; 

Zo= E 9u’'iQ)z^z'', (4.8) 

e(wi).(/i-j /)=0 


and Pfiu = R^^rnio)- We assume furthermore that Zq and Zi are real valued for / G L^(IR^,R^), 
and hence = g,y^ and = -G„fj,. 

With an appropriate canonical change of coordinates (that is, it preserves the r.h.s. of (4.2)) 
the term R_i and all non resonant terms in Rq and Ri cancel out. Indeed we have the following 
fact, which we quote from Theorem 6.4 [12]. 
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Proposition 4.3 (Birkhoff normal forms). There is a canonical transformation (z,/) —>■ 
where 

z' = z + n°’\Q,z,Q{f)J), 

f ^ ^ Si’i(Q,z,Q(/),/)). 

such that in the new coordinates (z^f) we have 

> 2,0 


(4.9) 


K{Q, z, /) = Qif)) + H^ + Ro + Ri+n with 

4 

^ Q(/), /) + E 0, /), 

i=2 


(4.10) 


with H 2 and Rj like in Lemma and where we have: 

(1) the term R^i in (4.4) is here R-i = 0; 

(2) all the nonzero terms in Rq with |^ + i/| < 2N + 1 are in normal form, that is e ■ {fi — v) = 0, 
and are in Zq; 

(3) all the nonzero terms in Ri with + < 2N are in normal form, that is e- [^ — v) £ ae{Rpg), 

and are in Z\; 

(4) we have R 2 £ C"*(U, C) for U C M x C” x K x PcZ-k 0 , neighborhood of [q(uji), 0, 0,0) and 


\R 2 {Q,z,f, e)| < C{\Q - q{oo^)\ + |z| + ||/||e_J||/||L,. 

Sketch of the proof. Our canonical transformation will be generated by Hamiltonian functions of the 
following form: 


E Ap,,.iQ,g)z>^z’^+ E z^^rn{BpAQ,Q)J)- 






The Hamilton vector flow generated from this Hamiltonian vector field will be 




( 


\fi+u\=m 


E '^jAp.AQ^O)^- + E Vj^-Q.{BpAQiQ)A) 

Zj Zj 


/(s)^/ + s E 




Thus, 


n n 

-Y^sAjAV)? + ‘^-^AB.J{Af{i)) - - Y^sAAzA + ‘2-A{c^j,f) 

i=i j=i 

+ E e-A-v)ApAQ^e)z^z''+ E e-A-V) BpAQ^Q)^)- 

\^Ay\—Tn \fiAy\—Tn—l 

e{li-iy)^0 e-(/i-;y)^(Te(.C„i) 


By the above, we can erase the nonresonant terms. 


□ 
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In (4.10) the functional K is written in a form which is essentially the same of (1.16). What 
follows in Sections 5-6 is rather close to the classical discussion in [7, 70]. The difference between 
these papers and [13, 3] lies in the fact that the latter two use in an essential form the Hamiltonian 
structure of the system to study higher order interactions between discrete and continuous modes. 
The ideas originate from [14]. The method [7, 70] does not work well with higher order interactions 
and requires very stringent restrictions on the spectrum of the linearization, which are completely 
eased in [13, 3]. 

As we have seen in the analysis of (1-16) the 2nd power, the structure of the Fermi golden rule 
will be easily seen in the framework provided by Proposition 4.3. However the informal analysis on 
(1.16) which we made in Section 1 has to be supplemented by a number of estimates, especially for 
the variable /. So we will need to write the equation for / and derive some estimates. 


5 Equations 

In the new coordinates {z,f) in Proposition 4.3 our NTS continues to be of the form (4.3). In 
particular we have 

/ = JV/7^2;°^(Q, Qif)) + z, Q{f), /) + JVfZoiQ, z, Q(f)) 

+ JVfZiiQ, z, Qif), /) + JVfTZiQ, z, Q(/), /), 

where 

V/(7^fc;l(Q, Qif)) + H!, + Zo) = C^J + A'Jf, 

A'■.= dQ^f)nl^iQ,Qif))+ Y. dQ^f)a^,iQ,Qif))z^T, (5.2) 

|/i+j/|>2 

e-(/i —i/)—0 

and similarly we split the 2nd line of (5.1) into 

A"J/-i ^ z^z'^G^,.iQ,Qif)) + JVfniQ,z,Q,f)\^^Q(f), 

(5-3) 

A" := [Zi(g, 2 , Qif),f) + TZiQ, z, Q(/), /)]. 

So finally we write the equation of / as 

f = Cui,^f + AJf - i Y, z'^z''G^^iQ,0)+ D\ where (5.4) 

A = A' + A" and 91 = V/7^(g, z, q, /)|,=q(/) - i (G^.(Q, Q(/)) - G^^iQ, 0)). 

We write the equations for z as 

= Jr(i?2 + ^o)+i i^,^(AG^.(g,g(/)),/) +(5.5) 

We set, for n the number associated to wi in (Hll), 

A = (Ai,...,An), (5.6) 

where the Xj = Aj(a;i) are introduced in (Hll). 
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Proposition 5.1. For cq sufficiently small there exists aC > s.t., given a solution u(t) of the NLS 
which satisfies supt>Qinf^gR l|w(t) — < e < co; for t £ I = [0,cxd) we have 


\\f\\L^ii,wff^) + \\f\\LfiiM'~‘) - admissible pairs {p,q), 

II' 2 ^IIl 2 ( 7 ) ^ Ce for all multi indices pL with X ■ fi > oji, 

||w - wiIIljo(j-) < Ce. 

Furthermore 


lim zit) = 0. 

t—¥-\-OC 


(5.7) 

(5.8) 

(5.9) 
(5.10) 


(5.11) 


Notice that the case {p,q) = (oo,2) in (5.7) and inequalities (5.9) and (5.10) are an easy 
consequence of supopinf^gR ||u(t) — < e. 

Proposition 5.1 implies Proposition 2.5. For the proof see [16], in particular Section 12. First 
of all, by standard arguments there is a /+ G iJ^(R^, C^) such that for the / in (5.7) and for the A' 
in (5.2) we have 

lim |l/(t)-e‘^(‘"^*+/o^'(^)'^^)e-‘^*^/+||j,i =0, (5.12) 

t^ + OO ^ ' 

see Lemma 7.2 [16]. We can express the solution u{t) as 


u{t) = e + P{oj{t)){z'j{t)f,j + z'jiffjfj + f'{t)) where 

= ^ + <;L(Q,^,0(/),/) and /' = e-'<.™(0>^-W)./)(/ + si;;„(Q, 2 ,Q(/),/)). 


(5.13) 


where {z, f) are the variables in Proposition 5.1 and (fc, m) are arbitrary. This follows from the fact 
that composing the change of variables (3.18) and (4.9) yields a change of variables like in (4.9), see 
[ 12 ]. 

(5.11) and (5.12) imply lim = 0 in 77^ and lim 77?’^ = 0 in R. 

t—>-+oo ’ t—^+00 ’ 

It is easy to see that when we plug (5.13) in ft = JVE(u) we get 

/ = J(-A + ^-7t°’2)/ + Gi(u), 


with Gi(it) € C°(Hf, L],). On the other hand, / satisfies also (5.1), which is of the form 

f = J (—A + wi + A')f + G2(m), 


with G 2 (u) G C°{F[f., LI.). Then we have 

K{u)f = Gi{u) — G 2 {u) with k{u) := wi — i) + 77.°’^ + A' . (5.14) 

We have k G G‘^(Ilf,R) and we claim that tt = 0. Indeed, if K(u(to)) 0 for a given solution, we 
can find solutions for which u„(t, •) G 5(K^), M„(to,-) u{tQ,-) in 77^(R^), ||Mn(to)||Li(R 2 ) —oo, 

Gj{un{to)) —>■ Gjiulfo)) and n{un{t[f)) —>■ K(u(to))- This yields a contradiction because on one 
hand ||u„(to)||Li(R 2 ) —>■ oo implies for the corresponding / coordinates ||/n(^o)llLi(R2) oo, on the 
other hand (5.14) is telling us ||/„(7 o)IIli(rO ^ \K{u{to))\-'^\\Gi{u{to)) - G 2 (u(fo))||Li(R 2 ) and so 

ll/n(7o)|lLi(R2) -ff OO. 
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Integrating wi — i? + 7^°’^ + A' = 0 and by lim 7^?’^ = 0 we get for a fixed i?o S K 

lim {ujit+ [ A'{s)ds —'d{t)) =-do. 

t-s-+oo Jq 

Then (5.12) for := e'^'^°f+ becomes 

lim ||e-^’’(‘)/(t)-e-^‘^/i+||^i=0. 

Using this in formula (5.13) we obtain 

lim ||u(t) - = 0 

t—>- + oo 


which yields (2.12). Combining (3.13) and (5.13) wee have 

q{uj{t)) + Q{f) + = q{uJo). 


Then 


lim q{uj{t)) = q{LOo) - Q{h+). 


Then (H5) implies that there must be a uj+ s.t. lim uj{t) = a;+. The last sentence of Proposition 

t^ + OO 

2.5 follows from (5.13) and the estimates in Proposition 5.1. 

□ 

By a standard continuity argument, Proposition 5.1 is a consequence of the following proposi¬ 
tion. 


Proposition 5.2. There exists eo > and a constant cq > 0 such that if T > 0 and if u{t) is a 
solution of the NTS which satisfies sup^gjinf., 5 gR l|w(t) — < e < eo where I = [0,T] then, 

if the inequalities (5.7)-(5.8) hold for this I and for C = Cq > cq, they hold also C = CqI2. 


6 Proof of Proposition 5.2 


The proof of Proposition 5.2 is basically the same of Proposition 6.7 in [16]. We give a schematic 
description of the main steps. The first is the following, which follows from theory in [21] and whose 
proof we review in Appendix A. 


Lemma 6.1. Assume the hypotheses of Prop. 5.2. Then there is a fixed c and an sq such that for 
all admissible pairs (p, q) and all s > Sq 


X-fL^LOi 

where we sum only on multiindices such that X ■ fj. — Xj < uji for any j such that for the j-th 
component of p. we have fXj 7 ^ 0 . 

□ 

The notation is simpler if we change frame. For M defined below we have M~^\JM = a^, see 
(1.17), and so we have 


:= M-HCu^M = aai-A + V + lj)+ 


where M := 


1 1 
—i i 


M-i = 


1 


2 11 -i 


0-3 = 
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( 6 . 2 ) 


Then we set ft. = M ^/, which satisfies for := M 0) and E ;= M 

ift = ft + Acr^h + ^ ( z^z^G^i, + E. 

e-(n-v)^a^[C^^) 


with A defined in (5.3). The last summation in (6.1) originates from the terms in (6.2) (or 

the corresponding ones in the 1st line of (5.4)). We cancel these terms by a normal forms argument 
For 

g = h + Y, Y:= ^ z'^z''R+^^iX ■ (u - g))G^,, (6.3) 

we have 


ig — ICi^^g + Aa^g + [iE — K,^-^Y\ + ^ ^ z^z'^G^i, + Aa^Y + E. 




We then compute 


iE = ^ X- (v — g)z^z''{X- {v — g))G^i, + T where 

\\-{n-v)\>u>i 

T ;= ^ [d^.Yigzj + XjZj) + - XjZj)] . 


(6.4) 


(6.5) 


Inserting (6.5) in (6.4) we obtain 

ig = K-uj^g + Aa^g + Aa^Y + T + E. 


( 6 . 6 ) 


So we have canceled the terms. Notice that T contains terms of this type but by (5.5) 

they are smaller. So g is smaller than ft. In fact the following is true, and is proved in Lemma 4.6 
in [21] (the statement in Lemma 4.6 [21] has a systematic typo and should be replaced by 

, where M there is like our s here). 

Lemma 6.2. Assume the hypotheses of Prop. 5.2. Then for fixed s > 1 there exist a fixed c such 
that if So is sufficiently small, for any preassigned and large L > 1 we have ||<?|| 2 , 2 ((q ^ ce. 

□ 

For the transpose of M, and using = 2M~^ , / = Mh and G^j, := M~^G^j,v{Q,t)), by 
direct computation we have 

i( JG^.(Q, 0), /) = i{M^JMM-^G^,{Q, 0), ft) = 2(fTia3G^,, ft). 

Notice that Gp,iy(Q,0) = —Gvfj.{Q,0), see in Lemma 4.1, implies by M ^ = aiM~^ 

G,^ = M"^G,^(g, 0) = -aiM-^G^^iQ, 0) = -aiG^,. 

Then substituting (6.3) in (5.5) we obtain 

iSj^j = ^jiH2 Y Zq) + 2 izj _ {g,cF^Gv^) + dzjR 

_ (6.T) 

+ 2 ^ V, -^-(i?+^^(A-(/3-a))G„^,a3G,^), 
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Let us consider the set of multi-indexes 

Ai := {a : \ ■ a > oji and X ■ a — Xk < V/c s.t. ak 0}. 
Set also A := {A • a : a S M}. 

Like in [16], there is a new set of variables C, = z + 0{z^) s.t. for a fixed C 

IIC-^IU? <CCoe2, IIC-^IUoo and 

s,iO = h) + a^.Zo(C, h) + V, 

+ 2 '^sGyo); 

\a=\v ^3 

where for a fixed constant cq we have 

n 

II^jCjIUmo.t] < co(l + C'o)e^. 

1=1 

Now we consider, like in [16], 

n n 

i=i i=i 

-4 ^ A-z/lm(^C“C''(^K,.,i(“^’“)^aO,cr3 G,.o 


A-Q: = A-J 


In the second line of ( 6 . 12 ) we have a sum 

r(C) := -4 E ilm (i?E 0, ^3 G(I:o) , for 

lga 

G(L,C):= E C“GaO. 


A-Q!=L 


( 6 . 8 ) 


(6.9) 


( 6 . 10 ) 


( 6 . 11 ) 


( 6 . 12 ) 


(6.13) 


For W = limj^oo e A+t^i)^ there exist _F(a,C) g 1 F^'P(]R 2 ^ ( 5 ^ 2 ^ Jqj. all fc g R andp € (l,oo) 

with 2G(L,C) = VFF(^’«), see [ 21 ]. Then for 

r(C) = -4 E 1™ Ini[(i?_A(-A - wi + ie)F^^’^\F['^’^^) - {Ra{u3i -L + ie)F^^’^\F['^’^^)] 

LeA^'^° 

= 4 ^ lim lm(i?A(wi -L + (6.14) 


lga 


= —4 lim 


LeA 

Now we assume: 


eNO 7r 2 (a:2 - (A - U3i))‘^ + e 




(H13) for some fixed constant F > 0 and for all ( g C" we have 
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(6.15) 


r(c) < -r ^ irr 

aeM 

Then integrating and exploiting (6.9) we get for t G [0,T] 

r lk“lli2(o,t) < cCoe^ - ^SjAj|zj(0)p+ ^SjAj|zj(t)|2. (6.16) 

a as in (H14) j j 

We want to conclude for some other fixed c' 

l.h.s. of (6.16) < Sc'Coe^. (6.17) 

Since supo<(<g^ inf^gR < e by hypothesis, we can conclude that A^ |zj(t)p < c'e^ 

for any t. Then (6.16) implies (6.17) since here we can assume Cq > 1 we get (6.17). 

We conclude that for cq > 0 sufficiently small and any T > 0, (5.7)-(5.8) in / = [0,r] and with 
C = Cq implies (5.7)-(5.8) in / = [0,r] with C = c(l + \/Cq) for c. This yields Proposition 5.2. 

□ 

Remark 6.3. Notice that in the proof of asymptotic stability in [13], where sj = — 1 for all j, orbital 
stability is a consequence of (6.17) rather than the other way around. Indeed in that case we have 

^ Aj|2;j(t)p + r ^2 ll^“llL2(o,t) ^ ^ Aj|2:j(0)|^, (6.18) 

j a as in (H14) j 

and taking the initial datum uq very close to we can assume - Aj|zj(0)p < c'e^. Then each 
term in the l.h.s. of (6.18) is small for all t > 0. Furthermore this and and the fact that from (5.5) 
we derive that the time derivatives Zj{t) remain small, we conclude that Zj{t) q fQj. 


7 Theorem 1.4 and cubic quintic equations 

We have not carried out numerical experiments to check examples to which Theorem 1.4 applies. 
In this section we combine numerical results in [63] with number of assumptions to propose some 
possible applications of Theorem 1.4. We discuss mainly hypothesis (H14). 

We consider the cubic-quintic NFS (1.5). For each m > 0 one can find a family of vortices 
ioi uj G O = (0,w*) with tpui > 0 [36]. Here, the upper bound w* of O is given by 
w* := supjo; > 0 | 3s > 0, + gs® < 0} which is in this case Notice that if 

^s^ — ^s^ + is® > 0 for all s > 0, by Pohozaev identity we can show that there are no nontrivial 
bound states in the energy space. 

As we have already discussed in Section 1, for each m = 1,2,3,4, 5 [63] shows that there is a 
critical value ujcr such that for uj > uicr the vortices are spectrally stable, while for cj < cocr they 
are spectrally unstable. For to = 1 the value is oJcr ~ 0.1487, see Section 5.5 [63]. In all these cases 
spectral instability is generated as follows. As lu approaches from above, two distinct eigenvalues 
on the imaginary axis iA^^^(a;) and iA^^^(w) coalesce at w = oJcr at a point iAcr (iAcr ~ i0.0478 for 
TO = 1). As w decreases further, two eigenvalues bifurcate from iAcr out of the imaginary axis. In 
[63] it is not stated explicitly whether or not the eigenvalues iA^t)[^;[) for j = 1,2 are simple and 
whether their algebraic and geometric dimensions coincide. The fact that only eigenvalues with 
different signatures can generate by collapse eigenvalues outside iK is well known, and we formalize 
it as follows. 

In the following, Co; will denote the linearized operator of the vortices of the 

cubic-quintic NFS (1.5). 
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Lemma 7.1. Consider equation (1-5), an m = 1,2, 3,4, 5 , a corresponding vortex and the operators 
Coj ■ Suppose that there exists Eq > 0 s.t. for w € (wcr, Wcr + Eo) eigenvalues (w) for j = 1,2 
are simple and their algebraic and geometric dimensions coincide (and so are 1). Then one of them 
satisfies (H 14 -). 

Proof. Recall, preliminarily, that if 2 ; is an eigenvalue of with z ^ —z, that is if z ^ iK., then 
il{v,v) = 0 for any v G kei(£uj — z)- This follows from 

z{J~^v,v) = {J~^CuiV,v) = -{v, J-^£ujv) = -z{J~^v,v), (7.1) 

which uses J~^Cui = If z = iA € R we have 

iX{J~^v,v) = iA(J“^u,U), 

so that {J~^£ujV,v) G R for any v G ker(£tj — iA). When iA £ R is simple we have {J~^£t^v,v) ^ 0 
for any v G ker(£^ — iA)\{0} and the sign is called the Krein signature. Since 

{J~^£^v,v) = \X{J~^v,v), 

it is clear that the Krein signature of iA is positive exactly if {J~^v,v) = i^ with <; < 0 (this 
explains why in (H14) the hypothesis Sj = 1 implies negative Krein signature and that positive 
Krein signature corresponds to Sj = —1). Notice that for iA to have a well defined Krein signature 
it is not necessary that be simple, and it is sufficient that {J~^£uiV,v) have constant sign as r; 7 ^ 0 
varies in ker(£(.j — iA). 

Now let us proceed with the proof of Lemma 7.1. By hypothesis the iA^'^^(u;) for w > ujcr are 
simple and hence they have a well defined Krein signature, which is constant in ui > ujcr (see [45]) and 

can be written in the form s^^'> := —ifl(^w \ ■C^^) for appropriate generators ^ G ker(£^ — iA('^^(w)). 
We have G { — 1,1}. For w < Wcr in [63] it is proved that there are two eigenvalues z^^^(a;) and 
zF)(aj) with Rez^'^^(aj) 7 ^ 0, which for reasons of symmetry satisfy z^^^(w) = —^*-^^( 0 ;). These two 
eigenvalues exit from iA*--’^ (tUcr) as w decreases. With an argument by contradiction we suppose now 
that fiF) = 3 ( 2 ) for Then, by Section 6.1 in [45], this continues to be true for w < uJcr in the 

sense that the quadratic form (J“^£,^-,”), which is definite in ker(£tj — iA^^)(a;))©ker(£(^ — iA(^)(a;)) 
for Lo > (jJcr must continue to be definite also in ker(£,^ — z^^^(a;)) © ker(£tj — z*'^^(a;)) for lo < oJcr- 
But from what we saw in (7.1) this is not possible. This gets us to a contradiction. So 7 ^ 
for UJ > Wcr and one of the two must be equal to 1 . □ 

Since (1.5) is translation invariant it is beyond the scope of our theory. In order to find 
examples of equations not translation invariant which satisfy (Hid) it is natural to add to (1.5) a 
small potential EK(|a;|) with a point of relative minimum in 0. We will show that the perturbed 
equation has vortices. As e —>■ 0 they converge to vortices of (1.5) in any space Efe(R^,C). 

Then the spectrum and the eigenfunctions of the linearizations £^'^ converge to spectrum and 
eigenfunctions of £uj. In particular, assuming our mixed rigorous and numerical proof that (1.5) 
satisfies (Hid) for w > Wcr, then we will have obtained this result also for the operators with 
E 7 ^ 0. 

Let (j)ui{C^r) = e^'^^ijjt^{r) be the vortex of [63] with > 0. Under the assumption that for 
(1.5) the kernel of £uj restricted in := {u G L'^ \ is radially symmetric} is Span} (in 

[63] it is shown that numerically this appears to be generically true) we show that for small e there 
exists = e™®V’i.cj,e(E) which is a solution of 

0 = + UJ(j)oj^e + £V4>bj,e ~ ( 7 - 2 ) 

More precisely, we have the following proposition. 
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Proposition 7.2. Assume ker£tjo |^2 = {J(j)uio}- Then there exist Jg > 0 and eg > 0 s.t. for 
£ € (—Eg,erg) and uj € (ajg — 5o,wo + dg) there exists 4>u),e G nfc>oSfc(R^,C) which satisfies (7.2). 
Furthermore, the map (w,e) —>■ 4>a;,s is in C^({u!o — do, tag + (^o) x (“£Oj£^o)) Sfc) for arbitrary k >0. 

To prove Proposition 7.2 we consider a preparatory and standard lemma. 

Lemma 7.3. Assume ker Cuj \]^2 = {J4>bj}- Then Au, = —Ar + ^ + uj — 3tjjf + is invertible in 
and < \\u\\l^. 

Proof. Let v € L^^^(R^,R) satisfy A,^v = 0. Then, multiplying by e‘"*® we have 
(—A + w — 2|()),^p + '3\(fui\^)u + {—(fl, + 2(^^|(^(^p) u = 0, 
where u = Therefore, using the natural identification between C and we have 

Cu,u = 0 . 

By the assumption, we have u = aJ(j)uj and thus v = aJipuj- However, v has values in R. x {0} while 
Jipuj has values in {0} x R. So a = 0 and kerH^ = {0}. Therefore, is invertible. Finally suppose 
A,^v = u. Then, first we have ||e‘"‘®u||i 2 = ||u||i 2 < ||u||i 2 . Next, multiplying by e*™®, we have 

(-A + u;)e'’"^ = {SiPl - 

Taking the norm of both sides, we have the conclusion. □ 

Proof of Proposition 1.2. Set d = oj — ivg and consider 4>ujg+s,e = {'f’uio+s + vs,e), where vs,e is 
radially symmetric and real valued. Then, substituting this into (7.2) and for v = vg.e, we have 

24^0^ = G{d,£,v) where G{d,£,v) := -£V{'fg,g+s + v) + u - dv + Ng,g+s{v), 

(7.3) 

with Ng,g.,gs(y) nonlinear in v and the convention [/(t)]|( = fib) — f{a). We can rephrase (7.3) as 

Fid,£,v) = 0 where Fid,£,v) := v — A~^Gid,£,v). (7.4) 

The function F is in C^((wg — i5g,wg + i5g) x (—Eg, eg) x i7^,iJ^). An elementary application of 
the implicit function theorem yields a function vs^e in C^((wg — (ig,a;g + (5g) x (—eg. Eg), iJ^). By a 
standard bootstrapping argument, FI^ can be replaced by for arbitrary k. □ 

In [63] it is checked numerically that for a generic vortex of the (1.5) 

AdgiTu)) — (jy^j, dx\4^uj, , d)u)(fu), J , J ^24^u)'\ • (^A) 

We have chosen H(|a;|) with a relative minimum at 0. The following fact is well known, see Theorem 
4.1 [24] and Theorem 3.0.2 [25]. 

Lemma 7.4. Consider an m = 1 vortex of the{1.5). Assume ^ 0 and (7.5) for the 

linearized operator Then for e > 0 sufficiently small we have 

ker£^^i = 5pa?r{ and Ng{C)fi) = Span{J<j)u,,e,dui(l)u:,e}- (7-6) 

Furthermore, has an eigenvalue of algebraic and geometric multiplicity 2 which is of the form 
i^(e) = iE-\/2e + o(e), with e the eigenvalue of the Hessian matrix of the potential in 0, and with 
eigenfunctions 

= y/2dxj4>u: + iVeeJxjcfuj + o(e), (7.7) 

with l|o(E)j|Ej, = o(e) for any k. 
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Proof. Everything is proved in Theorem 3.0.2 [25]. We only remark that there is only one small 
eigenvalue i/r(e) S iM+ which has to be of multiplicity 2. Indeed, implies by 

the symmetry J~^Vaj{xi,X 2 )J = Vuj{x 2 ,—x\) also J'^^i\x 2 t—xi) = i/r(£:)J'I'^'^^(x 2 , —xi). By 
J 4 >ui{xi,X 2 ) = - 4 >uj{x 2 ,-xi) we have 

J'^[^\x 2 ,-Xi) = jV 2 {da:.,(l)a,){x 2 ,-xi) - i^/eex 24 >u{x 2 ,-xi) + o(e) = 

V 2 dx 2 [jf>u>ix 2 ,-Xi)] +i^eJx24>ujixi,X2) + o(e) = V2dx24>ojixi,X2) + i^/eeJx24>u>{xl,X2) + o(e) 


and hence necessarily J'i’^^\x 2 , —Xi) = 'I' 2 ^^(a;i,a; 2 ). This implies that i^(e) has multiplicity 2. 


□ 


In the following remarks we discuss whether the perturbations of equation (1.5) satisfy the 
hypotheses of Theorem 1.4. 

Remark 7.5. (H1)-(H3) are trivially satisfied. 

Remark 7.6. It is known that equation (1.5) satisfies (H4) and the same holds for the perturbations. 
Remark 7.7. According to the numerical experiments in [63] hypothesis (H5) is true generically. The 
same will be true for the perturbations. 


Remark 7.8. (H6) is proved numerically for uj > LOcr in [63] and for the perturbations is a consequence 
of Proposition 7.2 and Lemma 7.4. 


Remark 7.9. Assuming the numerical results in [63] which claim that (7.5) is true for generic w, 
then for the perturbations (H7) is a consequence of Lemma 7.4. 

Remark 7.10. We don’t know the status of hypotheses (H8)-(H12) for (1.5), but if they hold, they 
hold also for the perturbations. 

As we have already mentioned, failure of (H8) would yield some Jordan block of dimension 2 or higher 
forcing 6*"^“ to grow algebraically in the invariant space A’^(£*) in (2.10) producing essentially a 
linear instability, which would yield an easy to detect nonlinear instability. Since this is not what 
the numerical experiments show, we conclude that probably [63, 77] confirm (H8). 

In [63] embedded eigenvalues and singularities at the edges are not discussed explicitly but hypotheses 
(H9)-(H10) seem to be confirmed. In [63] isolated eigenvalues are obtained as zeros of an Evans 
function. It is observed, see the discussion on pp. 371-372, that sometimes the Evans function 
is small near the continuous spectrum. This smallness is attributed not to eigenvalues sitting in 
the continuous spectrum, but rather to resonances on the other side of the continuous spectrum (in 
essence, to zeros of an analytic continuation of the Evans function beyond the continuous spectrum). 
Hypotheses (H11)-(H12) state that the eigenvalues between 0 and iw are positioned in a generic way. 
This is plausible to expect and probably can be proved for V generic. We did not attempt the proof. 

Remark 7.11. We have discussed at length Hypothesis (H13) in Section 1. It ought to be checked 
directly for perturbations of (1.5) or for equation (1.5). Notice that, since (1.5) is translation 
invariant, the search of an effective Hamiltonian is somewhat more involved, see [12, 16]. 

Remark 7.12. We have already discussed (H14) which is true, assuming that each of the eigenvalues 
(w) has some fixed Krein signature. Indeed this is what happens if each eigenvalue is simple, as 
we assumed in Lemma 7.1. Even if there is an eigenvalue of higher multiplicity, it is enough to ask 
for ifl{v,v) to have fixed sign for any eigenvector. This property on the signatures continue to hold 
for the perturbations. 
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A Appendix: proof of Lemma 6.1 

First of all, it is equivalent to consider equation 6.2 for h. Let Xc = and by an abuse of 

notation let us set Pc = M~^Pc{uji)M, where Pc is introduced under Lemma 2.3. Set also K. = JC^u 
The following three lemmas are Lemma 3.1-3.3 in [21]. 

Lemma A.l (Strichartz estimate). There exists a positive number C such that for any k € [0,2]; 
(a) for any h = Pch and any admissible all pair {p,q), 

(h) for any g(t,x) € and any couple of admissible pairs (jii,qi) (j> 2 ,q 2 ) we have 

II f e-'^*-^'>’^Pcg{s,-)ds\\^p, Kn < C\\g\\ 

Jo * ^ 

Lemma A.2. Let s > 1. 3 C = C such that: 

(a) for any f € S'(R^), 

< ClI/lU.; 

(b) for any g(t,x) G 


■)dt 


< C\\g\\r,rl 


LI 


Lemma A.3. Let s > 1. 3 C such that V g{t,x) G /S'(IR.^) and t G R; 


f e-'d^-^^’^Pcg{s,-)ds 

Jo 


L^Ll 


<C\\gl 


Lemma A.4. Let {p,q) be an admissible pair and let s > 1. 3 a constant C > 0 such that V 
g{tjx) G S'(R^) and t G R; 


Jo 


LfLl 




The following is Proposition 1.2 in [21]. 

Lemma A.5. The following limits are well defined isomorphism, inverse of each other: 

Wu = lim for any u G L^ 

t—>- + oo 

Zu = lim ^ 

t —^“|“00 

For any p G (1, oo) and any k the restrictions of W and Z to fl extend into operators such 

that for for a constant C we have 

ll^llw'“’P(R2),Wc'“’^ + W^Wwc ’^^ ^ 

with Wc'^ the closure in IF^’^’(R^) ofW'^'P{Mf) fl PcLl- 
The following is Lemma 3.5 [21]. 


31 














Lemma A.6. Consider the diagonal matrices = diagil^Q) E- = diag{0,l). Set P± = ZE±W 
with Z and W the wave operators associated to K,. Then we have for u = PcU 


1 

P+u = lim lim / + ie) — i?^(A — ie)] wdA 

e—)- 0 + 27 ri M^+oo J^ 

P-U = lim lim 

€—>■0+ 27ri M— >-+CX) 


/ — (jJ 

[i?K;(A + ie) - Ric{X- ie)] udX 

-M 


and for any Si and S 2 and for C = C(si, S 2 ) we have 


||(P+-P_-P,c73)/L2.., <C\\fh2.,,. 


Now we look at the term E in (6.2). 

Lemma A.7. For any preassigned s and for cq > 0 small enough we have 

E = i?i + i ?2 with ||i?i||Lj([o,T],ffi) + < C(s,C'o)e^. 

Furthermore for a fixed constant c we have 

II^I|l“((o,t),r) < cCqC^. 


(A.l) 


(A.2) 


Proof. The estimate on A = A' + A" follows from the definitions of A' in (5.2) and of A" in (5.3). 
E is a sum of various terms. For example we have 

||^^rM-i[G^.(Q, 0 )-G^.(Q,Q(/))]||^.([o_^],^i,.) 

< II^'^^1Il?[o,t]I|G^.( 0,0) - G^.(Q,0(/))IL^o.([o,^],^u«) < Cle\ 

So this term can be absorbed in R 2 . Another example is /3(|/P)/ = X|/|<i/3(|/P)/ + X|/|>i/3(|/P)/- 
The 1st term can be bounded, schematically, by 


lll/P/llL;([0.T].ffi) 


< 




I LB 


< 


LllO,T] - "LB{[0,T],Wf<>) 


< C^e^ 


while the 2nd term can be bounded by 


\w. 


1 . 2 I.||/||^ 2 L 


< 


\L-1 


< C^e^, 


wi 


(A.3) 


(A.4) 


where in the last step we use 


r r.-i 


< 


" Lir 


for some 0<a<lbyL>3 


(which we can always assume), interpolation and Sobolev embedding. 

Notice that by V/7^fe;^(Q, p, f)\e=QU) = Qif)^ f) have by (3.16) 

l|vX;L(Q,e,/)i.=Q(/)IL2ffi.^ < IIII/IL^,.-^IIl?(II/IIl= + kl + IQ(/)I)IIl? 

2^2 


- lu IIL?L 


2 r 2,-0 

t^x 


I^IUr + ll/IUr^')< 2 C^oe^- 

Consider for example the contribution of 


Vf [ BL{xJix),Q,z,gJ)f^ix)dxig=Q(^f)'^BL{xJ{x),Q,z,Q{f)J)f^ ^(x) 
+ [ deBLix, /(x), Q, z, Q(/), f)f^{x)dx + d 2 BL{x, /(x), Q, z, Q(/), f)f^{x). 


(A.5) 
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The last term can be treated like above, since \\Bl{x, f{x), Q, z, Q{f), /)||l“ < C by (3.20). We 
can use (A.3) or (A.4) for the 1st term of the r.h.s., since L — 1 > 3. Finally let us consider the 2nd 
term in the r.h.side. If we take g G L^H~^ we need to bound 

[ dt [ \{g,deBL{x,f{x),Q,z,Q{f),f))i^2 f^{x)dx\ 

Jo 

< [ dt\\\{g,deBL{x,f{x),Q,z,Q{f),f))L2^\\L2\\f^\\Li 

Jo 

< II 1(5, deBL{x, /(a:), Q, z, Q(/), /))i 2 ^ WJ^WlIli 
and we bound the last factor by (A.4). We have for fixed t 

|||(g, deBL{x, f{x),Q, z, Q{f),f))Ll, Wli < ||56Bl||b(e_,,eoI|5||s_, 

so that by (3.20), or by its analogue for the B^ in Lemma 4.1, we have that the last quantity is 
bounded by C||g||//-i. This yields a bound ||lst term 2nd line (A. 5 ) 11 ^ 17^2 < C^e^. 

Proof of Lemma 6.1 We rewrite (6.2) as 

iA=[/C/i + A(P+-P_)]/i + A[P,a3-P+ + P_]a3/i+ ^ + P^E. 

e.(/j-!^)Gcre(£^i) 


Then we have 


/i(t) =W(t, 0 )P*'^/i( 0 ) + / A[PcfT 3 -E++P_]CT 3 /r 

Jo L 


^ + PeE 


ds. 


e.(/i-l')eo-e(£ui) 

where the following operator commutes with /C: 

U{t,s) =^!lMs')d8fp+-p.) ^ 

Then 




< 


ll^(0)llffi + X/ ll^^^’^lli?II^Mi"llL“ 


Hi, 


(A.6) 


+ ll"4|Ut”ll^llL?ff,i-"" ll-^lllG-f^I Il'^2|li2j^i.«- 

The terms on the second line are O(e^) and the r.h.s. is bounded by the r.h.s. of (6.1), proving 
Lemma 6.1. 
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